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Abstract — Batteries have seen an increased use in electric ground 
and air vehicles for commercial, military, and space applications 
as the primary energy source. An important aspect of using 
batteries in such contexts is battery health monitoring. Batteries 
must be carefully monitored such that the battery health can be 
determined, and end of discharge and end of usable life events 
may be accurately predicted. For planetary rovers, battery 
health estimation and prediction is critical to mission planning 
and decision-making. We develop a model-based approach uti- 
lizing computaitonally efficient and accurate electrochemistry 
models of batteries. An unscented Kalman filter yields state 
estimates, which are then used to predict the future behavior of 
the batteries and, specifically, end of discharge. The prediction 
algorithm accounts for possible future power demands on the 
rover batteries in order to provide meaningful results and an 
accurate representation of prediction uncertainty. The frame- 
work is demonstrated on a set of lithium-ion batteries powering 
a rover at NASA Ames Research Center using real experimental 
field test data. 
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1. Introduction 

Batteries have seen an increased use as the primary energy 
source in electric ground and air vehicles for commercial, 
military, and space applications. In such contexts, battery 
health monitoring (BHM) becomes a critical issue, since 
failure of the battery directly affects the availability of the 
asset. Batteries must be carefully monitored such that the 
battery health can be determined, and end of discharge and 
end of usable life events may be accurately predicted. In 
planetary rovers, a battery health management (BHM) system 
is especially important, due to its role in short- and long-term 
mission planning and decision-making [1,2]. 

A BHM system consists of both battery state estimation 
and prediction components. While much research has been 
carried out on battery modeling and battery state estimation, 
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the prediction component, i.e., prognostics, has only recently 
begun to receive attention [3-5]. Even more recently, battery 
prognostics has begun to be applied to electric vehicles. 
In [6, 7], end of discharge prediction was developed for an 
unmanned aerial vehicle platform. In [8], an approach for 
predicting the remaining driving range of a battery-powered 
car was developed. These approaches differ by the type of 
battery model used, and the underlying estimation and pre- 
diction algorithms. Empirical battery models are used in [6], 
and equivalent circuit models are used in [7-9]. Equivalent 
circuit models are popular because they are relatively simple 
and computationally efficient. For estimation, the particle 
filter algorithm is used in [2, 3, 6], whereas the unscented 
Kalman filter (UKF) is used in [7-9]. The UKF is generally 
preferred over the particle filter when the model assumptions 
required for the UKF can be met, due to the much higher 
computational requirements of the particle filter [10]. Predic- 
tion algorithms used are typically sample-based algorithms, 
such as Monte Carlo sampling. They differ due to how 
future inputs are characterized and how that uncertainty is 
incorporated into the predictions [11]. 

Our approach to a BHM system is different from previ- 
ous approaches in several ways. First, we use a new 
electrochemistry-based battery model recently developed 
in [12]. Unlike similar models, it is composed only of 
ordinary differential equations, thus leading to a model that is 
as computationally efficient as equivalent circuit models. We 
use the UKF for state estimation, based on this new model. 
For prediction, we use an operator-centric approach that 
provides best-, average-, and worst-case usage predictions. 
The prognosis framework applied here is based on previous 
work presented in [9, 11], 

We apply our BHM system to a planetary rover testbed 
developed at NASA Ames Research Center, in which we 
monitor a set of 24 lithium-ion batteries that are used to 
power the rover motors [13]. We provide experimental results 
validating the rover battery model, and demonstrate battery 
state estimation and end-of-discharge prediction using real 
field test data. 

The paper is organized as follows. Section 2 describes the 
model-based prognostics framework. Section 3 summarizes 
the battery model. Sections 4 and 5 discuss the estimation 
and prediction approaches, respectively. Section 6 presents 
the experimental results. Section 7 concludes the paper. 


2. Model-based Prognostics 

In this section, we first formulate the prognostics problem, 
and present the uncertainty representation framework [11]. 
We then provide an architecture for model-based prognostics 
that will be applied for BHM. 
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Figure 1. Prognostics architecture. 


Problem Formulation 

We assume the system model may be generally defined as 

x(fc + 1) = f (k,x(k),0(k), u(fc), v(fc)), (1) 

y (k) = h(fc, x(fc), 0(k),u(k), n(fc)), (2) 

where fc is the discrete time variable, x(fc) £ R" x is the 
state vector, 0(k) £ is the unknown parameter vector, 
u(fc) £ is the input vector, v(fc) £ is the process 
noise vector, f is the state equation, y(/c) £ K" -1 ' is the output 
vector, n (k) £ is the measurement noise vector, and h 
is the output equation. 2 The unknown parameter vector 0(k) 
is used to capture explicit model parameters whose values are 
unknown and time-varying stochastically. 

Prognostics is concerned with predicting the occurrence of 
some event E that is defined with respect to the states, 
parameters, and inputs of the system. We define the event 
as the earliest instant that some event threshold T e : R" x x 
l" 8 x R n “ — y B, where B = {0, 1} changes from the value 
0 to 1. That is, the time of the event k E at some time of 
prediction kp is defined as 

k E (kp) = infjfc £ N : k > kp A T E (x.(k), 0(k ), u (k)) = 1}. 

(3) 

The time remaining until that event, A k E , is defined as 

A k E (k P ) = k E {k P ) - k P . (4) 

For batteries, the event E corresponds to end of discharge 
(EOD). EOD is defined by a voltage threshold V E od , where 
T e is defined by V < V E od- When the battery voltage is 
less than the cutoff voltage, V E od , EOD is reached and T E 
evaluates to 1. 

The prognostics problem is inherently uncertain, due to the 
random nature of the system evolution (due to v(fc)), and 
nondeterministic future inputs (u(fc) for k > kp). Therefore, 
k E and A k E are random variables, and we must compute 
p(k E (k P )\y(k 0 :kp)) and/or p(Ak E (k P )\y(k 0 :k P )) [9,14], 

Uncertainty Representation 

In order to predict k E , four sources of uncertainty must be 
dealt with in general [11]: (i) the initial state at time kp, 
x(fc„); (ii) the parameter values 0(k) for all k > kp, denoted 
as ®k P (the subscript kp indicates the start time of the 
trajectory); (iii) the inputs u(fc) for all k > kp, denoted 
as Ufc p ; and (iv) the process noise v(fc) for all k > kp, 
denoted as V /,, p . In order to make a prediction that accounts 
for this uncertainty, we require the probability distributions 
p(x), p(@k P ), p(Ufc P ), and p(V kp ). 


2 Bold typeface denotes vectors, and n a denotes the length of a vector a. 


For describing the probability distribution of a generic trajec- 
tory A*;, we introduce a set of surrogate random variables 
A a = [A 2 A 2 . . .] [11], We describe a trajectory using A a and 
instead define p( A a ), which in turn defines p( A k ). These 
surrogate variables can be used to describe trajectories in 
myriad ways. For the parameter, input, and process noise 
trajectories we have the surrogate variables A g, A„ , and A, . 
Additional discussion on the use of surrogate variables can be 
found in [11]. 

Prognostics Architecture 

We adopt a model-based prognostics architecture [11, 15], 
in which there are two sequential problems, (i) the esti- 
mation problem, which requires determining a joint state- 
parameter estimate p(yt(k),0(k)\y(ky.kp)) based on the 
history of observations up to time k, y(ko:kp), and (ii) 
the prediction problem, which determines at kp, using 
p(x(k),0(k)\y(k o :k P )), p(X e ), p(\ u ), andp(A„), a proba- 
bility distribution p(k E (kp)\y(ko:kp)). The distribution for 
A k E can be trivially computed ftom p(k E (kp)\y(ko:kp)) by 
subtracting kp from k E (kp). 

The prognostics architecture is shown in Fig. 1 [11]. In 
discrete time k, the system is provided with inputs rq r and 
provides measured outputs y^. The estimation module uses 
this information, along with the system model, to compute an 
estimate p(x.(k), 0(k)\y(ko:k)). The prediction module uses 
the joint state-parameter distribution and the system model, 
along with the distributions for the surrogate variables, p{ Xg), 
p(\ u ), and p(\ v ), to compute the probability distribution 
p(k E (kp)\y[ko:kp)) at given prediction times kp. 


3. Battery Modeling 

In order to apply model-based prognostics, we require a 
model. For this purpose, we employ an electrochemistry- 
based lithium ion battery model developed in [12]. In contrast 
to the empirical and equivalent circuit battery models used in 
previous approaches [3,8,9], the new model presented in [12] 
is based on the underlying electrochemical equations, but at a 
level of abstraction high enough that the model is still efficient 
with the improved fidelity. The model is represented as a 
set of ordinary differential equations and can be converted 
to a discrete-time representation and solved efficiently with a 
sample time of 1 s. We summarize here the model equations 
and refer the reader to [12] for additional details. 

The battery model computes the voltage as a function of 
time given the current drawn from the battery. Several 
electrochemical processes contribute to the cell’s potential. 
The different potentials are summarized in Fig. 2 (adapted 
from [16]). The overall battery voltage V(t) is the differ- 
ence between the potential at the positive current collector, 
<^s(0, £), and the negative current collector, <fi s (L,t), minus 
resistance losses at the current collectors (not shown in the 
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Table 1. Battery Model Parameters 





Figure 2. Battery voltages. 


diagram). 

The potentials at the current collectors are described by 
several voltage terms. At the positive current collector is the 
equilibrium potential Vu tP - This voltage is then reduced by 
V StP , due to the solid-phase ohmic resistance, and V„ p , the 
surface overpotential. The electrolyte ohmic resistance then 
causes another drop V e . At the negative electrode, there is 
a drop V v<n due to the surface overpotential, and a drop V s>n 
due to the solid-phase resistance. The voltage drops again due 
to the equilibrium potential at the negative current collector 
VJy n . We describe each of these voltages in the following 
subsections. 


Equilibrium Potential 

The equilibrium potential (also known as the open-circuit 
voltage) is captured using the Nernst equation: 

Vjj,i = Uq H — In ( + VlNT,i ; (5) 

nb \ Xi ) 

where i refers to the electrode (n for negative or p for 
positive), Uq is a reference potential, R is the universal gas 
constant, T is the electrode temperature (in K), n is the 
number of electrons transferred in the reaction (n = 1 for Li- 
ion), F is Faraday’s constant, x is the mole fraction of lithium 
ions in the lithium-intercalated host material [17]. During 
discharge, Li ions move out of the negative electrode and into 
the positive electrode, so x n decreases while x p increases, 
and Vijjj — Vjj, n will decrease. 


Here, Vint,* is the activity correction term (0 in the ideal 
condition). We use the Redlich-Kister expansion: 


Vint,* — — -p 
nb 



(Ot — 11 fc+1 — Xi) 

1 * j (2^-1)!-^ 


Here we use N p = 12 and N n = 0. The identified parameters 
are given in Table 1 [12]. 


We let qi represent the amount of Li ions in electrode i, as 
measured in Coulombs. The flow of Li ions is opposite to the 
flow of current, so (/,; changes in the same direction as x t , and 
we have: 


- Ch 
Xi gmm ’ 


(7) 


Parameter 

Value 

U 0 ,p 

4.03 V 

A p , o 

-33642.23 J/mol 

Ap,i 

0.11 J/mol 

Ap,2 

23506.89 J/mol 

Ap,3 

-74679.26 J/mol 

Ap, 4 

14359.34 J/mol 

Ap, 5 

307849.79 J/mol 

Ap, 6 

85053.13 J/mol 

Ap,7 

-1075148.06 J/mol 

Ap,8 

2173.62 J/mol 

Ap,9 

991586.68 J/mol 

Ap, io 

283423.47 J/mol 

Ap,n 

-163020.34 J/mol 

Ap, 12 

-470297.35 J/mol 

Uo,n 

0.01 V 

-^-71,0 

86.19 J/mol 


where r/ max = q p + q n refers to the total amount of available 
Li ions. It follows that x p + x n = 1. When fully charged, 
x p = 0.4 and x n = 0.6. When fully discharged, x p = 1 and 

x n = 0. 3 

Concentration Overpotential 

When a battery is discharged, the reactions take place at the 
surface of the electrode and this results in a concentration 
gradient across the cell. The model accommodates this by 
splitting the total electrode volume into two individual control 
volumes (CVs), one for the bulk (with subscript b) and one for 
the surface (with subscript s). 

For the volumes, the concentration of Li ions is computed as 


Qb,i 

,z — 

(8) 

Vb,i 


Qs,i 

Cs.i j 

(9) 

V s ,i 



where, for CV v in electrode i, c V: i is the concentration and 
v Vt i is the volume. We define Vi = v b,i + v s ^. Note now that 
the following relations hold: 


Q p — Qs, P ~b Qb, P ( 10 ) 

Qn = Qs,n Qb,n ( 11 ) 

Q = Qs, P ~b Qb, P ~b Qs,n ~b qb : n- ( 12 ) 


As the battery discharges, Li ions move from the surface layer 
at the negative electrode, through the bulk, and to the surface 
layer at the positive electrode, in order to match the flow of 
electrons. Li ions also move from the bulk CV to the surface 
due to the concentration gradient. The diffusion rate from the 
bulk to the surface is expressed as 

Qbs,i — (Cb,i (13) 

where D is the diffusion constant. 


3 For LLC 0 O 2 , x must be at least 0.4; Li cannot be reversibly removed 
beyond that [17], 
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The q variables are described as 


Qs, p I'app T Qbs,p 

Qb,p Qbs,p “b 'Lapp ^app 

Qb,n Qbs,n "F t app t app 

Qs,n t app “b qbs,m 


(14) 

(15) 

(16) 
(17) 


where i app is the applied electric current. Initially, c, vi = c/ v , 
thus diffusion is zero. 


The concentration overpotential is the difference in voltage 
between the two CVs due to the difference in concentra- 
tion. Using the expression for equilibrium potential, we can 
compute the potential for the bulk volume and the potential 
for the surface layer; the difference between them is the 
concentration overpotential. So, we can explicitly account 
for the concentration overpotential simply by using as the 
expression for equilibrium potential, the equilibrium potential 
of the surface layer [12], i.e.. 


Vu,i = Uq 


RT 

nF 


In 



+ VlNT.i, 


(18) 


State of Charge 

The state of charge (SOC) of a battery is defined to be 1 when 
the battery is fully charged and 0 when the battery is fully 
discharged by convention. In this model, it is analogous to 
the mole fraction x n , but scaled from 0 to 1. We distinguish 
here between nominal SOC and apparent SOC [12]. Nominal 
SOC is computed based on the combination of the bulk and 
surface layer CVs in the negative electrode, whereas apparent 
SOC is be computed based only on the surface layer. When 
a battery reaches the voltage cutoff, apparent SOC is 0, and 
nominal SOC is greater than 0 (how much greater depends 
on the difference between the diffusion rate and the current 
drawn). Once the concentration gradient settles out, the 
surface layer will be partially replenished and apparent SOC 
will rise while nominal SOC remains the same. Nominal (n) 
and apparent (a) SOC are defined using 


SOC n 

SOC a 


Qn 

0.6 g max 

Qs,n 

0.6 ’ 


(25) 

(26) 


where g maXs ." = 9 max ^ 4 


where x s ^ is computed using 


Qs,i 


(19) 


and 


max max 

y s.i y 


Vi 


( 20 ) 


Ohmic Overpotential 

The voltage drops due to the solid-phase ohmic resistances, 
the electrolyte ohmic resistance, and the resistances at the 
current collectors are constant and lumped together into re- 
sistance R 0 : 


l'o — lappRo- 


( 21 ) 


Surface Overpotential 

The overpotentials due to charge transfer resistance and SEI 
kinetics are described by the Butler- Volmer equation, which, 
for Li ion batteries, reduces to 

K,, = ^arcsinh ( ' JOj , (22) 

where Ji is the current density, and ./,q is the exchange 
current density. The current densities are defined as 

J* = ^ (23) 

Ji 0 = ki(l - x s>i ) 0 ‘{x Sti ) 1 ~ a , (24) 


Battery Voltage 

Battery voltage can now be expressed as follows: 

v = V U}P -V Uin - V' - V’ >p - v' tn , (27) 

where 

v' = {Vo - Vf)/T 0 (28) 

Kp = ( V n,P - <p)/tv,p ( 29 ) 

Kn = (Y v , n - K, JKn, (30) 

and the r parameters are empirical time constants (used since 
the voltages do not change instantaneously). 

The model contains as states x, q SiP , q b}P , q b<n , q s ^ n , V', V p p , 

and V^ n . The single model output is V. Identified model 
parameters are given in Table 2 [12], 

Model Validation 

Measured and predicted voltage curves for a 0.044 A dis- 
charge (approximately equal to open-circuit voltage) are 
shown in Fig. 3, and for a 2 A discharge in Fig. 4. The data 
here was obtained from a laboratory setting using the rover 
batteries. Clearly, the model is very accurate. Measured and 
predicted voltage curves for battery data from a rover field 
test are shown in Fig. 5. The model is still quite accurate in 
this case. Note that the spike in predicted voltage occurring 
near 4000 s is due to a gap in the recorded data, causing an 
incorrect sample time to be given to the model. 


4. Estimation 

In order to accurately predict the future behavior of the 
battery, we must first estimate its state. For this purpose, 
we use the unscented Kalman filter (UKF) [18, 19]. Among 
nonlinear filters, the UKF generally has better accuracy than 


where fcj is a lumped parameter of several constants including 
a rate coefficient, electrolyte concentration, and maximum 
ion concentration. 


4 Note that SOC of 1 corresponds to the point where q n = 0.6 , since 
the mole fraction at the positive electrode cannot go below 0.4, as described 
earlier. 
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Table 2. Battery Model Parameters 


Parameter 

Value 

Q maX 

1.32 x 10 4 C 

R 

8.314 J/mol/K 

T 

292 K 

F 

96487 C/mol 

n 

1 

D 

7.0 x 10 6 mol s/C/m 3 

To 

10 s 

a 

0.5 

Ro 

0.085 n 

S p 

2 x 10 -4 m 2 

kp 

2 x 10 4 A/m 2 

V s,p 

2 x 10" 6 m 3 

Vb,p 

2 x 10" 5 m 3 

T V,P 

90s 

S n 

2 x 10" 4 m 2 

kn 

2 x 10 4 A/m 2 

^Sj-n 

2 x 10" 6 m 3 

V b,n 

2 x 10" 5 m 3 

Tp^n 

90s 



Figure 3. Comparison of predicted and measured open- 
circuit voltage. 



Figure 5. Comparison of predicted and measured voltage for 
a rover battery during a field test. 


random variable y by some nonlinear function y = g(x), 
and computes the mean y and covariance P yy using a set 
of deterministically selected weighted samples, called sigma 
points [18]. X‘ denotes the v'th sigma point from x and 
w 1 denotes its weight. The sigma points are always chosen 
such that the mean and covariance match those of the original 
distribution, x and P xx . Each sigma point is passed through 
g to obtain new sigma points y, i.e., 


y* = g(Af‘) 

(31) 

with mean and covariance 


y = Y^w i y 

(32) 

Vyv = Y. wi{ y i -y)(y i -y) T - 

(33) 


I 


We use the symmetric unscented transform, which selects 
2n x + 1 sigma points symmetrically about the mean [19]: 



Figure 4. Comparison of predicted and measured voltage for 
a constant 2 A discharge. 


the extended Kalman filter, and avoids the high computational 
cost of particle filters [20]. We summarize the filter basics 
here; more details may be found in [18, 19]. 

The UKF approximates a distribution using the unscented 
transform (UT). The UT takes a random variable x £ R rix , 
with mean x and covariance P xx , which is related to a second 


X 1 


K 

(n x + k) 1 
2 {n x + k) ’ 


i = 0 

* = 2 n x 



(34) 


(35) 


where y\/(n x + k)P xx J refers to the ith column of the 

matrix square root of {n x + n)P xx (e.g., computed using the 
Cholesky decomposition). The number k is a free parameter 
that can be used to tune the higher order moments of the 
distribution, and if x is assumed Gaussian, then selecting 
k = 3 — n x is recommended [18]. 

The UKF assumes the general nonlinear form of the state 
and output equations, but is restricted to additive Gaussian 
noise. First, n s sigma points X k _i\ k _ 1 are derived from the 
current mean Xfc_ 1 ifc_ 1 and covariance estimates Pfc_ 1 i^_ 1 . 


5 


The prediction step is: 


^ k\k— 1 l|fc— 1 > U k— 1) 5 ^ l,...,77- s (36) 

fc|fc— 1 l)? ^ !?•••? (37) 

n s 

* k \k-i = J2 wiX Uk-i ( 38 ) 

i 

n s 

y klk -i = j2 wi yUk-i ( 39 ) 

i 


Pk|k-1 — Q+ 

n s 

E ^(*k|k-i - *fc|*-i)(^’fc|fc- 1 - * fc | fe _i) T , (40) 

i 

where Q is the process noise covariance matrix. 

The update step is: 

n s 

P yy = P + E ^CVfcik-i _ yfc|k-l)CVfc|k-i _ yfc|fc-l) T 

(41) 

n s 

p x V = E wI (^fcik-i - *k\k-i){yi\k-i - y k \ k ~i) T 

(42) 

Kfc = P^P" 1 (43) 

Xfc|k = x fc | fc _i +K fc (y fc - yfcife-r) (44) 

Pfc|k = Pfc|k-i - KfePyj/Kki (45) 

where R is the sensor noise covariance matrix. 

For the battery, we have 7 states, resulting in 15 sigma points. 
An estimate of SOC can be computed from the estimates of 
the states. This can be computed using the unscented trans- 
form, in which the sigma points for the state are transformed 
into sigma points for SOC, using the equations for computing 
SOC from the battery charge. We can then extract mean and 
variance of the SOC estimate from these transformed sigma 
points. 


5. Prediction 

Prediction is initiated at a given time kp using the current 
joint state-parameter estimate, p(x(kp), 0(kp)\y(ko'-kp)). 
The goal is to compute p(kp(kp)\y(k 0 :kp)) using the state- 
parameter estimates and assumptions about uncertainty re- 
garding the future parameter, input, and process noise values. 

In this work, we assume all parameters are known, so 0(k) 
is empty. Further, we assume that, because the model is 
so accurate, process noise is negligible. Of course, process 
noise is not zero, but in this case the uncertainty in the future 
inputs dominates significantly, so including process noise has 
virtually no effect on the prediction results. In the following, 
we describe the approach for the general case as originally 
developed in [11]. 

For one realization of each of the uncertain quantities at 
prediction time kp: the state x(fcp), the parameter trajectory 


Algorithm 1 k E (k P ) <- P(x(fc P ), Q fcp ,U fcp , Vfc P ) 

1: k <— kp 

2: x(fc) <— x(fcp) 

3: while Ts(x(fc),0 fcp (fc),U fep (/c)) = 0 do 

4: x(fc + 1) «- f (k, x(fc), © fcp (k), U fcp (fc), V fcp (fc)) 

5: k «— k + 1 

6: x(fc) <— x(fc + 1) 

7: end while 

8: fcp(fcp) <— k 


® kp , the input trajectory XJ kp , and the process noise trajec- 
tory V kp , the corresponding realization of kp can be com- 
puted with the system model as shown in Algorithm 1 [11]. 
In Algorithm 1, the function P simulates the system model 
until the threshold Tp evaluates to 1. 

This algorithm requires computing first realizations of the 
state-parameter distribution, the parameter trajectory, the in- 
put trajectory, and the process noise trajectory. The distribu- 
tion for the state comes from the UKF, and the distributions 
for the parameter, input, and process noise trajectories are 
defined indirectly by the set of surrogate variables. So, we 
are interested in computing the distribution for kp from the 
distributions for p(x(fc P ), 0(fc P )), p(Xg), p( X u ), and p( X v ) 
(see Fig. 1). 

In this paper, since only future input uncertainty is con- 
sidered, we need to define only the surrogate variables for 
the future input trajectory. Here, future input trajectories 
may take many complex forms. Instead of representing 
these complexities directly, we consider instead an equivalent 
constant-loading distribution for the future inputs. That is, 
we assume that the future power usage for a battery will be 
constant with the value drawn from some distribution. 

Prediction methods differ by how they sample from the 
given distributions, and how they call P. Different methods 
were investigated in [11], including Monte Carlo sampling, 
unscented transform sampling, and the inverse first-order- 
reliability method (FORM). In the case of the rover, the 
operator really only needs to know EOD predictions for best- 
, average-, and worst-case usage scenarios. For the state 
estimate, we use as samples the sigma points provided by the 
UKF. Each sample is simulated forward three times, once for 
each use case. From this we obtain best-, average-, and worst- 
case EOD predictions, each with some small variance (due to 
the state estimate variance). 

If instead the operator desires a probability distribution, then 
an equivalent constant-loading distribution can be selected. 
We can determine the appropriate statistics for an equivalent 
constant-loading distribution by analyzing field experiments. 
Assuming a normal distribution, we can compute the average 
power for each field experiment, then compute the mean and 
variance of the average power. For a particular realization 
of the input, we can sample the average power from this 
distribution, and use that power for k > kp. In this paper 
we use a mean power draw of 3.3 W with a variance of 1.0. 

In this case, we have one surrogate variable needed to de- 
scribe the future input, which is used to define the constant 
power draw, as defined by the determined statistics. In 
order to sample efficiently from this distribution, we use the 
unscented transform method, as originally described in [9]. 
The unscented transform described in Section 4 is used to 
sample from the future input distribution, and in this case 
only 3 samples are needed since there is only a single 
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Figure 6. Battery current. 



surrogate variable. This is much more efficient than the 
hundreds or thousands of samples required for Monte Carlo 
sampling. Prediction then proceeds in the normal way, with 
the P function being called for the sigma point combinations 
(each state sigma point is simulated for each future input 
sigma point). Using the unscented transform equations we 
can then determine the statistics of EOD, namely, the mean 
and variance. (For additional statistical moments, extended 
versions of the unscented transform or Monte Carlo sampling 
are required. Alternatively, the inverse FORM method can 
be used to determine cumulative distribution function values 
at desired points.) Mathematical details of this approach are 
provided in [9, 1 1]. 


6. Results 

In this section, we present experimental results obtained 
from the planetary rover testbed at NASA Ames Research 
Center [13]. The rover motors are powered by 24 lithium-ion 
batteries, with two strings of 12 batteries in series, connected 
in parallel. So, each battery sees only about half of the total 
current required to operate the motors. Fully charged, each 
battery supplies 4.2 V. The batteries can be safely discharged 
to 2.5 V, but the drop from about 3.3 V down to 2.5 V happens 
relatively quickly, so to prevent overdischarge, the voltage 
cutoff used is Veod = 3.3 V. 

The BHM system monitors the health state of each battery, 
and periodically provides EOD predictions, to let the rover 
operator know how much longer the rover can be used in the 
field. We demonstrate the BHM system as it performs on a 
field test of the rover. In this scenario, the rover was executing 
various maneuvers, e.g., moving straight, turning at various 
speeds, and stopping. The current drawn from the batteries 
during these maneuvers is shown in Fig. 6. 

As described in Section 3, the battery model is very accurate. 
As a result, the battery state can be confidently estimated, as 
shown in Fig. 5. The SOC estimate for a single battery is 
shown in Fig. 7. We plot here the apparent SOC, which is 
computed with respect to the charge available at the electrode 
surface. Due to the development of concentration gradients, 
apparent SOC may increase when the diffusion rate from the 
bulk to the surface layer exceeds the discharge rate of the 
battery, as observed in the figure. On larger time scales, SOC 
is seen to decrease steadily, and rover operation ceases once 
SOC reaches around 5%. 

The future power demands on the battery are unknown, espe- 
cially in the driving scenario presented here, where an opera- 
tor is deciding spontanteously where to drive the rover. If the 
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future usage of the rover is known exactly, the BHM system 
should be able to obtain accurate predictions, since the model 
is accurate. Prediction results in this case are shown in Fig. 8, 
and the results are very accurate. There is some slight positive 
bias, which is due in part to sensor noise, because EOD is 
determined based on measured battery voltage, which, due 
to noise, will be earlier than actual EOD because the voltage 
threshold will be reached first due to noise. Overall, results 
are still very accurate, and have virtually no spread, since the 
estimation results had very little variance. Relative accuracy 
of the predictions averages to 98%, with a spread of only 
0.16% relative standard deviation. Accuracy can only be 
improved with a better model, and, given the model we have, 
represents the best possible performance. This serves as a 
baseline for the scenarios where the future inputs are not 
known. 

If the future battery usage is unknown, we provide the op- 
erator with prediction results for three separate cases, best-, 
average-, and worst-case power demands. Prediction results 
in this case are shown in Fig. 9. Clearly, predictions assuming 
average-case usage are not accurate, as the actual power 
usage is less than the average case up until 4000 s; after that 
the average power usage of the remainder of the scenario is 
greater than the average case. The average relative accuracy 
in this case is 83.3% based on the predictions corresponding 
to the average case. The predictions based on average usage 
still falls within the bounds set by the best- and worst-case 
usage. These results are with respect to continued usage, so 
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Figure 9. Aks predictions at best-, average-, and worst-case 
usage. 
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Figure 10. Akp: predictions for equivalent constant-loading. 


the operator knows how much driving time is remaining, i.e., 
the predictions do not include time spent while the rover is 
stopped. 

Using the equivalent constant-loading distribution, we obtain 
the predictions given in Fig. 10. Here, because the actual 
mean power drawn in this scenario is larger than in the 
collection of field experiments, the median A kg predictions 
are larger than the true values. Relative accuracy computed 
using these values is only 63%. Because of the large variety 
of scenarios in the field experiments, the spread is large, 
at 45% relative standard deviation. These numbers can be 
improved only if additional information is known about the 
future usage and that information is captured in the future 
input characterization, e.g., through the use of additional 
surrogate variables [11]. 


7. Conclusions 

In this paper, we described a battery health monitoring 
framework for a rover. The prognostics framework was 


detailed, along with algorithms for state estimation and end- 
of-discharge prediction. The provided battery model was 
validated with both laboratory and field test data. Validation 
of the BHM system was demonstrated using data from an 
actual field test of the rover. 

Future work will include end-of-life prediction for the batter- 
ies, which must track the degradation of the battery health 
over multiple usage cycles. Initial results suggest that the 
battery model described here is applicable in this case [12]. 
In addition, when there is some knowledge available as to 
the future operation of the rover, e.g., a set of waypoints, 
or a large set of past driving scenarios, future inputs can be 
characterized in a more complex way than presented here, 
which can potentially lead to more accurate results with much 
less spread [11]. However, this depends on what kind of 
information is desired by the operator, as in many cases, 
predictions for best-, average-, and worst-case usage are 
sufficient. 
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